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Abstract. Wc calculate the limiting gap distribution for the fractional parts of log n, 
where n runs through all positive integers. By rescaling the sequence, the proof quickly 
reduces to an argument used by Barra and Gaspard in the context of level spacing 
statistics for quantum graphs. The key ingredient is Weyl equidistribution of irrational 
translations on multi-dimensional tori. Our results extend to logarithms with arbitrary 
base; we deduce explicit formulas when the base is transcendental or the rth root of 
an integer. If the base is close to one, the gap distribution is close to the exponential 
distribution. 

§1. The gap distribution is a popular measure to quantify the degree of randomness 
in a given deterministic sequence. Rudnick and Zaharescu have shown [T2] that the gap 
distribution of the fractional parts of 2 n a (n = 1, . . . , N), for almost all a, converges in 
the limit iV — > oo to an exponential distribution — the gap distribution of a Poisson point 
process. (2 n may be replaced by any lacunary sequence.) The same is expected to hold 
for the fractional parts of n k a, for any integer k > 2 and a of bounded type, but this 
remains unproven [101 [HI El H] • Numerical experiments also suggest that the fractional 
parts of n 13 may have an exponential gap distribution provided /3e(0,|)U(|,l). In the 
case /? — |, Elkies and McMullen [3] have calculated an explicit formula for the limiting 
gap distribution, which is evidently not exponential. Remarkably, their distribution co- 
incides with the limiting gap distribution of directions of points in the affine lattice Z 2 +a 
for any fixed vector a ^ Q 2 [9j. Both of these findings follow from the equidistribution 
of translates of certain curves on the homogeneous space ASL(2, Z)\ ASL(2, R). In the 
present paper we show that also the fractional parts of log b n have a non-exponential 
limiting gap distribution, which we calculate explicitly. Our derivation reduces quickly 
to an argument used by Barra and Gaspard [T] in the context of spectral statistics of 
quantum graphs. The key ingredient here is Weyl equidistribution on multi-dimensional 
tori, similar in spirit to [5J Chap. 3]. 

§2. To state our main results, let us denote by £ n the fractional parts of log b n. We 
denote by {£ n ,iv}n=i,...,iv the ordered set of the first N elements of £ n , so that 

(1) o < a,iv < 6,iv < ■ • • < Un < i. 

For purely notational reasons it will be convenient to set £n+i,n '■= £l,jv + 1 5 an d also 
consider the gap between the first and last element mod 1; this additional gap has of 
course no effect on the existence and form of the limiting gap distribution. 

For given N, the gap distribution Pn(s) is defined as the probability density on R>o, 

1 - 

(2) P N (S) = JjJ2 5 ( S - N (^+l,N - £n,iv)). 

n=l 
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Figure 1. The distribution of gaps between the fractional parts of logn, 
where n = 1, . . . , 10 4 . The piecewise continuous curve is the limit distri- 
bution P(s) of Theorem [I] for b = e. 




Figure 2. The distribution of gaps between the fractional parts of log b n 
with base b = e 1/5 = 1.221402758 . . ., where n = 1, . . . , 10 4 . The blue 
curve is the limit distribution P(s) of Theorem [TJ the red curve is the 
exponential distribution e~ s . 



We denote by T the set of transcendental numbers 6>l;6 = eisa natural example, 
cf. Figure [TJ The technique we present here also works for algebraic b, but in general 
leads to more intricate limit distributions. In ^10 we discuss the simple case of integer 
base b, and in §[Tl~1 the case when the base is the rth root of an integer. 



Theorem 1. Let b G T . For any bounded continuous f : M>o — > M, 

poo poo 

(3) lim / f(s)P N (s)ds= f(s)P(s)ds, 

N ->°° Jo Jo 
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where 

(^ b s- 2 (F(s\ogb) - F(sb-Hogb)) i fO<s<^- b 

(4) P(s)=l-^ b s-*F(sb-Hogb) if^<s<ik> 

(5) F(a) = a-^-(a; b' 1 ) - (a; b' 1 ) 

da 

and 

oc 

(6) (a; q) = J^[(l — aq n ) (q-Pochhammer symbol). 

n=0 

The derivative of the Pochhammer symbol is explicitly 

d °° i 

i=o y 

Note that P(s) has jump discontinuities at s = and however both the left and 

right limits exist at these points, and are finite. We also note that lim,,_ 5 . + P(s) = r^r ■ 
A more elementary formulation, which is equivalent to the above, is that for every 
s > 

(8) lim —\{ n <N: N{£ n+ltN - > s}\ = / P(s') ds' . 



We will show in the final paragraph Sjl2| of this paper that, in the limit b — > 1, P(s) 
converges to the exponential distribution e~ s , cf. Figure |2j 

§3. To gain a better insight into the limiting gap distribution P(s) for fixed b, let 
us recall that the fractional parts of log b n are not uniformly distributed mod 1 [6j. 
However, the sequence of r\ n (n = 1, . . . , N) given by the fractional parts of log 6 (n/iV) = 
log fe n — log fe iV has a limit density mod 1, and evidently the same gap distribution (to 
see this, recall that we have added the gap between the first and last element of the 
sequence mod 1, and the ordering in the sequence remains the same). It is an easy 
exercise to show that for any interval Ac. [0,1], 

(9) lim \z\{n < N : r] n G A}\ = [ p(x) dx 

JV-foo N J A 

with density 

(10) p{x)= l ^-b\ 



§4. A more refined statement of Theorem [T] is to consider the joint distribution of rj n 
and the subsequent gap; define the corresponding probability density on [0, 1] x M> by 

1 N 

(11) Pn(x,s) = j^^2S(x -r) njN ) 6(s-N(r] n+hN -r] n)N )), 

n=l 

where {^ n ,Ar}n=l,...,Ar is again the ordered set of the first N elements of r\ n . The joint 
distribution allows us to compare the statistics of gaps between elements near one point 
x of the unit interval with the gap statistics at a second point x' . 
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, i — 7: (a: b 1 ) if < a < 1 



The explicit formula for the limit distribution will involve the density 

(12) R(a, b) = (IT 1 ; b' 1 ) 8{a - 1) + I da 

[0 ifo>l. 

The second derivative of the Pochhammer symbol is 

(13) W M = (a;9) £,(^-a)(g-*-a)' 
Theorem 2. Lei 6 G T. For any bounded continuous f : [0, 1] 

/•I /*oo /*1 /*oo 

(14) lim / / f(x,s)P^(x,s)dsdx = / / f(x,s)P(x,s)dsdx, 
N ^°°Jo Jo Jo Jo 

where 

(15) P(x, s) = [b*' 1 log b) 2 R{sb x - x log 6, 6). 

§5. To unravel the previous statements further, let us consider a strictly increasing 
function [0, 1] — > [0, 1] which maps the sequence rjj to a new sequence fjj which now is 
uniformly distributed mod 1. This function is given by (cf. [7]) 

f x b x - 1 

(16) x^ J p(y) dy= - - - , 

and thus 

6^-1 - 1 

(!') Vn = -T — , ?7n,JV 



We have, explicitly, for n e [iV6~ fc , iV6~ fc+1 ) and k e Z> 



b - I 6 - I 

x ) and 
b k n — N 



(18) *• N (b-iy 

We will see in ^8] that the knowledge of the joint distribution Pn(x, s) for the original 
sequence and the analogue for the rescaled and ordered sequence, 

1 - 

(19) P/v(x, s) = — ^2 <K X _ Vu,n) S(s- N(fj n+1>N - fj n>N )) 

n=l 

are equivalent: 

Theorem 3. Let b G T ■ For any bounded continuous f : [0, 1] 

/>1 poo pi poo 

(20) lim / / f(x,s)PN(x,s)dsdx= / / f(x,s)P(s)dsdx, 
*-**>Jo </o Jo Jo 

(21) P(s) = (l-6- 1 ) 2 ^((l-6" 1 )s»6)- 

Note that the limit distribution is independent of x. This is not a consequence of fjj 
being uniformly distributed mod 1, but reflects the fact that the gap distribution is the 
same if we restrict the sequence to an arbitrary subinterval of [0, 1]. 

We will first prove Theorem [3] in £j6]and §[7J and then retrace our steps by subsequently 
showing in £j8] that Theorem [3] implies Theorem [2j where 

(22) P(x,s) = P{p(x)s)p{x) 2 = (b*- 1 log b) 2 R{sb x - x log 6,6). 
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We obtain Theorem [T] from Theorem [2] by taking an x-independent test function. 

§6. We now turn to the proof of Theorem[3]and study a more general setting. Consider 
a finite sequence .. . ,uij of positive real numbers (the "frequencies"). With each Uj 
we associate a real number /3j. We are interested in the statistical properties of the 
multiset 

j 

(23) S= l+IO^ + u^Z). 

i=i 

Here l±l denotes multiset sum, i.e. a union where we keep track of the multiplicity of 
each element (this is only relevant if the sets (3j + wJ 1 Z, j — 1, . . . , J are not mutually 
disjoint). The number of elements of S that fall into the interval / = [t — ^, t + ^ ] is 

j 

(24) N(t,L) = J2Nj(t,L), Nj(t,L) = |/D (/?,• +uJ 1 Z)\. 

3=1 

Evidently Nj(t,L) is periodic in £ with period a;" 1 . Furthermore 

(25) WjL - 1 < Nj(t, L) < ujL + 1. 
Therefore 

(26) hm^%^ = y^, 

3=1 

which means that the asymptotic density of «S in R equals ■ Uj. 

For any x G R we write |x| z G [0, |] for the distance between x and the nearest 
integer. We also set 

(27) n{x,L) = | [x - f,x+ f] nZ|. 
The verification of the following lemma is a simple exercise. 

Lemma 1. Given x G R, L > and G Z> 0; u>e /iawe n(x,L) = k if and only if 
\x — is both > \{k — L) and > |(L — k). 

In particular, if x is picked at random with respect to the uniform probability measure 
P on [0,1], then 



(28) P (n(z, L) = k) = E^k, L) :-- 



l-\k-L\ ifL-l<k<L + l 
otherwise. 



The following observation is due to Barra and Gaspard [1J . 

Theorem 4. Assume that the frequencies ui, . . . ,uj are linearly independent over 
Then, for — oo < a < b < oo and k G Z>o, 

v 7 t^oo (b-a)T K J 

where 

j 

(30) ^ n^^'^)- 

fei+...+fcj=fc j=i 
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In particular, 

n/=i t 1 - u j L ) i f L < max(wi, • • • , ujjY 1 



(31) E(0,L) 



otherwise. 



The convergence in (29) uniform over all choices of fix, . . . , f3j G E keeping 
fixed). 

Proof. Decomposing [aT, bT] into intervals of length u^ 1 and using N(t, L) = YLj=i Nj(t, L), 
we see that meas({t G [aT, bT] : iV(i, L) = k}) equals 

(32) 

n b -l 

^2 E ^ meas ({ s e [°» !) : ^(a/f^s + n),L) = kj for j = 1, . . . , J}) + E, 

k\+...+kj=k n=n a 

where n a = [u>\aT\, n b = [uibT\, and E is a real number whose absolute value 
is bounded above by the measure of the symmetric difference between [aT, bT] and 
[n a u^ 1 , n b uj^ ]; thus < 2u;^ 1 . Let us define a^ n G E/Z through the relation 
/3j + wJ^Z = nwjf 1 + coj (aj >n + Z). Then Nj(co^ (s + n),L) = kj holds if and only 
if Ifwf 1 ^ — ^u^s + f ] H uj l {otj^ n + Z)| = fcj, and by Lemma [l] this holds if and 
only if 2|^-s — o>j >n — is both > kj — UjL and > UjL — kj. It follows that 

meas({t G [aT, bT] : N(t, L) = k}) equals 

n b -l 

(33) Yl E % 

) + E, 

ki+...+kj=k n=n a 

where / Al ,..., fe>7 : (M/Z) J -> [0, 1] is defined by 

(34) f ku ... :kj (a>i,a 2 ,. . . ,a,j) = \\l[2 —s - otj - -kj > [kj - UjL\ ) ds. 

Jo ~~[ v w i z z / 

Note that each function fk ly ..,kj is continuous on (IR/Z) J . It follows from our definition 
of that a^.n+i = ctj ; n — co/wT/ 1 mod Z. In particular ai_ n is independent of n; in 
fact a\ t n = on '■= wiPi in R/Z for all n. On the other hand our assumption implies that 
the numbers l,U2^i i^z^x > • • • > ^j^i" 1 are linearly independent over Q, and hence by 
Weyl equidistribution we have 

/ oc \ Sn=% /fei,...,fcj( Q; l,n) ^2,71, • • • , Q!j,n) /" „ , N , , 

(35) ► / f ku ... :kj {a 1 ,a 2 ,...,aj)da 2 ...daj 

n b — n a JQBL/zy- 1 

as % - n a — > oo. Note that the convergence here is uniform over all choices of 
(3i, . . . , (3 j G E. The right hand side equals 



(se) / / n j ( 2 



— s — a,- fc,- 

wi J 2 J 



> — cjjLl ) (is, 



and changing order of integration and then substituting aj := Xj + ^-s — \kj (xj G 
IR/Z) we see that the expression factors as Ylj=i E\(kj,ujjL). Finally noticing also that 
n b — n a ~ Wi(b — a)T as T — > oo, we obtain (29). □ 

We now turn to the gap distribution. We order the elements of S and label them as 

(37) . . . < A_ 2 < A_i < A < Ai < A 2 < . . . . 
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Theorem 5. Assume that the frequencies Uj are linearly independent over Q, and that 
U\ is the largest among Ui, . . . ,uj . Then, for — oo < a < b < oo, and any s > 0, 
s ^ wf 1 , 

1 f 00 

(38) lim — \{j EZ: \j E [aT,bT\ and X j+1 - Xj > s}\ = / P w {s')ds' 

t— >oo (o — a)T ' J s 



where 



(39) P UJ ( s ) = uj 1 l[(l-^y(s-^ 1 ) + 



3=2 



< h,i=l 





n (!-^ s ) 



if < s < u 1 



if s > uj 1 1 . 



-i 



TTie convergence in (38) uniform over all choices of fii, . . . , f3j G M f&irf keeping 
U!,...,uj fixed). 

Proof. This is a corollary of Theorem |1J cf. e.g. Theorem 2.2 in [7j. The limiting gap 
density is given by the formula 



(40) 



d 2 E(0,s) 
ds 2 



□ 



§7. Proof of Theorem^ It suffices to prove that ( p0| holds when / is the characteristic 
function of a box [0,a] x [0,A] with A ^ (1 — b^^^ 1 , i.e. to prove that for any fixed 
a G [0, 1] and A > 0, A ^ (1 - Jr 1 ) -1 , we have 



(41) 

as — 7- oo, where 



G aA {N)^a / P(s)ds 



(42) G a>A (N) := i | {n G {!,..., iV} : ^ jJV g [0,a], TV^+i,^ - ^jv) G [0,A]} 



(cf., e.g., Example 2.3]). 
Set wj = 6-^(6-1) and /3 J (7V) = 



for j = 1, 2, . . .. For any given J, iV G Z>i 



we note that by (18), {Nfji, . . . ,Nrj N } equals the multiset sum of (f3^ +co~ 1 Z) n [0, N) 
for j = 1, 2, . . . , J and a remaining multiset of cardinality \Nb~ J ] consisting of Nfj n for 
n = N and all n < Nb~ J . Hence if we order the elements of Sn '■= W^ =1 (/3^ iV ' 1 + uJ lr L) 
as ... < A_2 < A_i < Ao < Ai < A2 < . . ., then 

(43) G a>A {N) = e Z : X j e [0,aN], X j+1 - A, g [0, + O^- 1 + r J ). 

Since 6 is transcendental, the frequencies U\, . . . , Uj are linearly independent over Q; 
thus Theorem [5] applies, and we obtain, for any fixed J G Z>i, 



(44) 



lim sup G ajA (N) < a / P w ,j{ s ') ds ' + °( 6 ~ 

N^oo Jo 

liminf G oA (iV) >a I P u ds' - 0(b~ 

AM-oo J Q 
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where P Wj j(s) = if s > u 1 1 while for s < u 1 1 we have 

(45) P UiJ (s) =vif[(l- %)8{8 - u?) + ^ II (! - <*j*) • 

i=2 1 i=i 

We now compute that 
(46) 



fo 



if A > u x 1 
WjA) if A < u^ 1 . 



7=1 i=l 



Hence, letting J — >■ oo in (44) and using Uj = b 3 (b — 1), we conclude 
(47) 







ifA>^ 



lim G atA (N) = a + a < , Q . ,.\ 

W_((i_ r i )s;r i))^ ifA <^ 



(4 / -P(s) cfe, 
'0 



where the last equality follows by a direct computation using (12) and (21). Hence (41) 
holds, and Theorem [3] is proved. □ 

§8. Let us now prove that Theorem [3] implies Theorem [2] (An entirely analogous 
argument also shows that Theorem [2] implies Theorem |3j so that the statements of 
these two theorems are in fact equivalent.) Given a bounded continuous function / : 
[0, 1] x ]R> — y K, we wish to prove that (14) holds. Since the sequence of probability 



measures Pn(x, s) ds dx is tight (cf., e.g., Lemma 2.1 in [7]), a familiar approximation 
argument shows that without loss of generality we may assume that / has compact 
support, i.e. there is some B > such that f(x, s) — whenever s > B. 

Let T : [0, 1] -> [0, 1] be the map in dEb and define T : [0, 1] x M> -> [0, 1] x M> 



through %(x,s) = (T(x),T'(x)s). This is a bijection, with inverse given by 

s 



(48) 



I -1 1 



x, s 



x , 



'T'CT-^a:))/' 



Now apply Theorem |3] with the test function / ol 1 . Then in (20) we have 

rl poo j 

(49) / / f o c Z~ 1 (x,s) P N (x,s) dsdx = jz^f(rj n , N ,s n> , 



J0 



n=l 



where s n . 



iv 



T(r) n+ x tN ) - T(r] n , N ) 



Let J-7v be the set of those n G {1, . . . , A^} for 



T '(Vn,N) 

which ?7 n+ i i Ar — r] n N < BbN^ 1 . Then since T'(x) = p(x) is continuous and bounded 
away from zero on [0, 1] we have s n ^N = N(r) n +i^ — Vn,N) + o(l) uniformly over all 
n G J-jv as N — > oo, and hence, since / is uniformly continuous, 

(50) f{Vn,N, s h ,n) = f(Vn,N, N(r] n+1>N - r] ntN )) + o(l), 

uniformly over all n G J-jy. On the other hand if n ^ J-jy then N(rj n+ltN — r] njN ) > 
Bb > B and also s n>y y > Bb = B, so that both sides of (50) vanish. Hence (50) in 



fact holds uniformly over all n G {1, . . . , iV} as iV — > oo, and it follows that the limit 
in the left hand side of (14) exists, and equals the limit in (20) (with / o T -1 ). Finally 
substituting (x, s) = %(x', s') in the right hand side of (20) we see (via (J22[)) that this 
limit equals the right hand side of (14), and the proof is complete. 
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Figure 3. The distribution of gaps between the fractional parts of log 6 n 
with base b = 10, where n = 1, . . . , 10 4 . The piecewise continuous curve 



is the limit distribution P(s) in (56). 



§9. Theorem [T] is now a direct consequence of Theorem [2] As to the explicit formula 
for the gap distribution, 

-i 



P(s) 



(51) 



P(x, s) dx 

[ {b x ~ l \ogb) 2 R{sb x - l \ogb,b)dx. 
Jo 



With the variable substitution a = sb x 1 log6, da = sb x 1 (log6) 2 dx = (logb)adx we 
have 

fs log b 



(52) 



Pis) 



s 2 log b 



R(a, b) a da. 



logb 



Recall the definition of R(a,b) in (12). The integral over the first term in (12) yields 
^logfe ^ < 1 < s log b and otherwise. For the second term, we use integra- 

tion by parts, 



(53) 



a— —-(a; b~ l ) da 

ao 9a2 



-i ai 



«i 



a^M- 1 ) 
d 

a— (a; b' 1 ) - (a;^ 1 ) 



d_ 
da 



-l Ol 



a; 6 ) da 



where ao = min{l, sb 1 logo} and a\ = min{l, slog 6}. This yields Q. 

§10. As we have noted, the technique presented here works also for algebraic 6. Let 
us consider briefly the simplest case of b > 1 being an integer. In this case, the multiset 
sum Sn = ^j = i(^ NS> + ujJ lr L) considered in the proof of Theorem [3] has all its points 
lying in (3^ + w^Z, i.e. the only gap lengths appearing are and o;^ 1 = (1 
furthermore Sn is periodic with period cuj 1 = b 



j-i 



LO-, 



One finds that the average 
multiplicity of (3^ + to^k in Sn for varying k E Z equals \_ b -t ■ Hence following the 
proof of Theorem [3] and taking the limit J — > 00, we conclude that for b > 1 an integer, 
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Figure 4. The distribution of gaps between the fractional parts of log 6 n 
with base b = a/To, where n — 1, . . . , 10 4 . The piecewise continuous curve 



is the limit distribution P(s) in (62) with r = 2 




Figure 5. The distribution of gaps between the fractional parts of 
101og 10 n (= log fe n with base b = 10 1/10 = 1.258925411...), where 
n = 1,...,10 4 . The blue curve is the limit distribution P(s) in (62) 
with r = 10. The red curve is the exponential distribution e~ s . 



we have a limit result as in (20) in Theorem |3| but with the limit density being given 
simply by 



(54) 



p(s) = b-^is) + (i - b- 1 )s(s - (i - r 1 )- 1 ). 



As to the limit of the unsealed sequence r] n , the first equality in (22) yields 

b l-x' 



(55) 



P(.r. .s) = bT-His) + (logft)^" 1 *! s - — 
o — l \ logo 
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and therefore 

f -i- s- 2 if -±- < s < jL - 
(56) P(s) = 6 _1 <J(s) + <^ logfc logb logb 

I otherwise. 

§11. The findings of the previous paragraph can be generalised to bases of the form 
b = m 1 ^, with m and r positive integers such that m l l p is not an integer for any prime 
divisor p of r. In this case, the multiset sum Sn = \±}j =1 (fi( N > + coJ 1 Z) considered in 
the proof of Theorem [3] has all points in the set U^ =1 (/3^ + uJ l X). Our assumption 
implies that the polynomial X r — m is irreducible over Q (since any monic polynomial 
of degree d < r dividing X r — m must have zeroth coefficient of absolute value m d//r , 
which is irrational). Hence Ui, . . . ,u r are linearly independent over Q, and we now see 



by the same argument as in the proof of Theorem |4j that the limit (29) exists in the 
case k = and is given by 



YY j=l (l-b-i(b-l)L) ifL<(l-6- 1 )- 1 
otherwise. 



(57) E(0,L)-- 
Using the finite q-Pochhammer symbol 

■r-l 

1 1 -aq n " 




if L < (1 - fe- 1 )- 1 
otherwise. 



We note that for L < (1 



r-l 



(60) -^-E(0 } L) = (l-b- 1 )E(0 1 L)J2 



j=0 y 



and thus 



(61) -l L E(,,L) 



r-l 

= < 1 -»- 1 )E^ = i 



The fact that this value is less than 1 is due to the non-trivial multiplicity of values 
in our sequence, which results in a positive density of zero gaps. The statements of 
Theorems [lfj3] therefore hold with the following limit distributions. The limiting gap 
distribution of Theorem Q] is 

f ^- b s-\F r {s\ogb) - Friab-Hogb)) if < a < ^ 

(62) P(s)=b- r 5(s)+l-^- b s- 2 F r (sb-Hogb) if ^- b <s<^- b 

with 

(63) F r (a) = a-^-(a;b- 1 ) r -(a;b- 1 ) r . 
The joint limiting gap distribution of Theorem [2] is 

(64) P(x, s) = b x - r 5{s) + (b x ~ l log bfRrisb*- 1 log b, b) 
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and the rescaled gap distribution of Theorem [3] is 

(65) P{s) = b~ r 5(s) + (1 - b^fRr ((1 - b~ x )8, b) 
with 

(66) R r (a, b) = (6 _1 ; 6 _1 ) r _i <J(a - 1) + < da 2 





a; 6 x ) r if < a < 1 
if a > 1. 



§12. We conclude this investigation with a few remarks on the family of limit distri- 
butions which we have obtained for different values of b G T. It is helpful to introduce 
the random point process in R given by the sequence 



(67) 



(b) 



{N(rj n + t + m) : 1 < n < N, m G Z} 



where t is uniformly distributed in [0,1]. This process is clearly stationary and has 
intensity one. By using the same strategy as in the proof of Theorem[3j one can deduce 
from TheoremU a limit law for the one- dimensional distribution of cp^ . For every closed 
interval A C R, k G Z>o, 



(68) 
where 

(69) 



lim ] 



\<p { 3 n A| 



fc) = E {b \k,me&s(A)), 



E (b) (k,L) = lim V TT^ife-,^- l)6 _j L) ; 



fcl+...+fcJ=fc j=l 



with Ei(k,L) as in ( |28[ ). In the case k = 0, relation (68) in fact follows directly from 
Theorem [3] (exploiting again Theorem 2.2 in [7]), where 



(70) 



pQ>), 



d 2 E^(0,s) 
di 2 



Although (68) requires b G T, the family of distributions (69) is well defined for any 



real b G (1, 00). Let us briefly analyse the limiting cases b — > 00 and b — > 1. 

For & ->■ 00, E^(k,L) converges to the statistics of the point process given by the 
randomly shifted integer lattice Z + t, with t uniformly distributed in [0, 1]. That is, 
E ( - b \k,L) — > E\{k,L) as in (28). The rescaled gap distribution satisfies 



(71) P {b \s) ^R(s, 00) = 5(s- 1). 

The limit of the raw p( b \s) distribution is slightly more involved since the asymptotic 
density p( b '(x) in (10) converges to 8(x — 1). We therefore cannot expect P^(s) to 
converge to a non-trivial limit without further rescaling. Let us extend to a 



probability density on 
(72) 



by setting p 



(b)( 



for x > 1. We have the scaling limit 



log b 



P 



(6) 



X 



log b 



-»■ e" 



which, however, is not a probability density on R>o- Working directly from Q, we see 
that for s fixed, 



(73) 



log b 



p(b) 



log b 



— >■ 







if < s < 1 
if s > 1. 
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In the limit 6 -»• 1, E^(k,L) converges to the Poisson distribution. That is, 

T k 

(74) E®(k, L) ->■ =-e" L , P {b) (s) ->■ e~ s . 
To see this, note that for L fixed and b sufficiently close to 1, 

oo 

(75) log E {b) (0, L) = M 1 - ( b ~ l)b~ j L) = -L + 0(b - 1) , 

i=i 

which proves the claim for = and the gap distribution P(s) (via Theorem 2.2 in [7]). 
Again, for L fixed and b sufficiently close to 1, we have for the remaining cases k > 1, 

(^ L )- E (°» L ) .2^ El(0 , (ft - 1)b -n L) . . . El (0, (b - l)b~i*L) 

(76) =(6-iw>(o,l)_ (i-^-i^r ^ 



Jl<...<J* 



(1 - (6- 1)&-J'*L) 



L fc EW(0,L) ^1 + 0(6-1)^) 



As to the limit of the unsealed sequence r] n , we have for the density 
(77) P {b \x)^l, x 6 [0,1], 



i.e., we have uniform distribution mod one. In view of d22b, PW(x, s) ->■ e" s and thus, 
for the raw gap distribution, P^)(s) — >■ e~ s . 

The above discussion can be readily adapted to bases of the form b = m 1//r considered 
in £11 In particular, we again observe that both the raw and rescaled limiting gap 
distributions (62) and (65) converge to the exponential distribution e~ s when r — > oo 
for fixed m, cf. Fig. [5| 

Acknowledgments. We thank Jon Keating and Zeev Rudnick for their comments on 
the first draft of this paper. 

References 

[1] F. Barra and P. Gaspard, On the level spacing distribution in quantum graphs. J. Statist. Phys. 
101 (2000) 283-319. 

[2] P. Billingsley, Convergence of probability measures, John Wiley & Sons Inc., New York, 1999. 
[3] N.D. Elkies and C.T. McMullen, Gaps in y/n mod 1 and ergodic theory. Duke Math. J. 123 (2004) 
95-139. 

[4] D.R. Heath-Brown, Pair correlation for fractional parts of an 2 . Math. Proc. Cambridge Philos. 
Soc. 148 (2010) 385-407. 

[5] M. Kac, Statistical independence in probability analysis and number theory. The Carus Mathe- 
matical Monographs, No. 12. Published by the Mathematical Association of America. Distributed 
by John Wiley and Sons Inc., New York 1959. 

[6] L. Kuipcrs and H. Niederreiter, Uniform distribution of sequences. Pure and Applied Mathematics. 
Wiley-Interscience, New York-London-Sydney, 1974. 

[7] J. Marklof, Distribution modulo one and Ratner's theorem. Equidistribution in number theory, 
an introduction, 217-244, NATO Sci. Ser. II Math. Phys. Chem., 237, Springer, Dordrecht, 2007. 

[8] J. Marklof and A. Strombergsson, Equidistribution of Kronecker sequences along closed horocycles. 
Geom. Funct. Anal. 13 (2003) 1239-1280. 

[9] J. Marklof and A. Strombergsson, The distribution of free path lengths in the periodic Lorentz 
gas and related lattice point problems. Ann. of Math. (2) 172 (2010) 1949-2033. 
[10] Z. Rudnick and P. Sarnak, The pair correlation function of fractional parts of polynomials, Comm. 
Math. Phys. 194 (1998) 61-70. 



14 



JENS MARKLOF AND ANDREAS STROMBERGSSON 



[11] Z. Rudnick, P. Sarnak and A. Zaharescu, The distribution of spacings between the fractional parts 

of n 2 a. Invent. Math. 145 (2001) 37-57. 
[12] Z. Rudnick and A. Zaharescu, The distribution of spacings between fractional parts of lacunary 

sequences. Forum Math. 14 (2002) 691-712. 

School of Mathematics, University of Bristol, Bristol BS8 1TW, U.K. 
j .marklof@bristol.ac.uk 

Department of Mathematics, Box 480, Uppsala University, SE-75106 Uppsala, Sweden 
astrombeOmath .uu.se 



